Once a reservoir, always a reservoir? Seasonality affects the pathogen maintenance potential of amphibian hosts.

Host species that can independently maintain a pathogen in a host community and contribute to infection in other species are important targets for disease management. However, the potential of host species to maintain a pathogen is not fixed over time, and an important challenge is understanding how within- and across-season variability in host maintenance potential affects pathogen persistence over longer time scales relevant for disease management (e.g., years). Here, we sought to understand the causes and consequences of seasonal infection dynamics in leopard frogs (Rana sphenocephala and Rana pipiens) infected with the fungal pathogen Batrachochytrium dendrobatidis (Bd). We addressed three questions broadly applicable to seasonal host-parasite systems. First, to what degree are observed seasonal patterns in infection driven by temperature-dependent infection processes compared to seasonal host demographic processes? Second, how does seasonal variation in maintenance potential affect long-term pathogen persistence in multi-host communities? Third, does high deterministic maintenance potential relate to the long-term stochastic persistence of pathogens in host populations with seasonal infection dynamics? To answer these questions, we used field data collected over 3 years on >1400 amphibians across four geographic locations, laboratory and mesocosm experiments, and a novel mathematical model. We found that the mechanisms that drive seasonal prevalence were different from those driving seasonal infection intensity. Seasonal variation in Bd prevalence was driven primarily by changes in host contact rates associated with breeding migrations to and from aquatic habitat. In contrast, seasonal changes in infection intensity were driven by temperature-induced changes in Bd growth rate. Using our model, we found that the maintenance potential of leopard frogs varied significantly throughout the year and that seasonal troughs in infection prevalence made it unlikely that leopard frogs were responsible for long-term Bd persistence in these seasonal amphibian communities, highlighting the importance of alternative pathogen reservoirs for Bd persistence. Our results have broad implications for management in seasonal host-pathogen systems, showing that seasonal changes in host and pathogen vital rates, rather than the depletion of susceptible hosts, can lead to troughs in pathogen prevalence and stochastic pathogen extirpation.

Evaluating environmental DNA as a tool for detecting an amphibian pathogen using an optimized extraction method.

Environmental DNA (eDNA) detection is a valuable conservation tool that can be used to identify and monitor imperiled or invasive species and wildlife pathogens. Batrachochytrium pathogens are of global conservation concern because they are a leading cause of amphibian decline. While eDNA techniques have been used to detect Batrachochytrium DNA in the environment, a systematic comparison of extraction methods across environmental samples is lacking. In this study, we first compared eDNA extraction methods and found that a soil extraction kit (Qiagen PowerSoil) was the most effective for detecting Batrachochytrium dendrobatidis in water samples. The PowerSoil extraction had a minimum detection level of 100 zoospores and had a two- to four-fold higher detection probability than other commonly used extraction methods (e.g., QIAamp extraction, DNeasy+Qiashredder extraction method, respectively). Next, we used this extraction method on field-collected water and sediment samples and were able to detect pathogen DNA in both. While field-collected water filters were equivalent to amphibian skin swab samples in detecting the presence of pathogen DNA, the seasonal patterns in pathogen quantity were different between skin swabs and water samples. Detection rate was lowest in sediment samples. We also found that detection probability increases with the volume of water filtered. Our results indicate that water filter eDNA samples can be accurate in detecting pathogen presence at the habitat scale but their utility for quantifying pathogen loads in the environment appears limited. We suggest that eDNA techniques be used for early warning detection to guide animal sampling efforts.

Optimized Batrachochytrium dendrobatidis DNA extraction of swab samples results in imperfect detection particularly when infection intensities are low.

Accurate detection of the amphibian fungal pathogen Batrachochytrium dendrobatidis (Bd) is critical for wildlife disease research; however, false negatives in detection do occur. Here we compared different DNA extraction methods to determine the threshold for Bd detection and identify an optimal extraction method to improve detection and quantification of the pathogen. We extracted both lab-created cell suspension standards using PrepMan Ultra, Chelex resin, and 3 spin column DNA extraction kits (Qiagen DNeasy Blood and Tissue, Zymo Quick DNA miniprep, and IBI gMAX mini kit), and further compared extraction methods using field-collected samples. We found that when extracting Bd DNA from cells in lab-created culture, the spin column extraction methods and PrepMan Ultra were equivalent, while the resin method detected higher Bd DNA quantities, especially at higher loads. However, when swabs from live animals were analyzed, low Bd quantities were more than twice as likely to be detected using a spin column extraction than with the PrepMan Ultra extraction method. All tested spin column extraction methods performed similarly across both field and lab samples. Samples containing low Bd quantities yielded inconsistent detection and quantification of Bd DNA copies regardless of extraction method. To manage imperfect detection of Bd, we suggest that presence/absence analyses are more informative than attempting to quantify Bd DNA when quantities are low. Overall, we recommend that a cost-benefit analysis of target species susceptibility and epidemiology be taken into consideration when designing an experiment to determine the most appropriate DNA extraction method to be used, because sometimes detecting low Bd quantities is imperative to the study, whereas in other situations, detecting low DNA quantities is less important.

Heterozygosity predicts clutch and egg size but not plasticity in a house sparrow population with no evidence of inbreeding.

We investigated the link between heterozygosity and the reaction norm attributes of reproductive performance in female house sparrows (Passer domesticus). We collected data on clutch size, egg size, hatching success and nestling survival in 2816 nesting attempts made by 791 marked individuals over a 16-year period. Pedigree analysis revealed no evidence of inbreeding. Neither parent-offspring regression nor an animal model revealed significant heritability in clutch or egg size. We selected 42 females that laid at least seven clutches at our study site and used a survey of 21 autosomal microsatellite loci to estimate heterozygosity for each female. We controlled for phenotypic plasticity and found that both clutch and egg size showed significant positive correlations with heterozygosity. We found no evidence that heterozygosity influenced the slope of individual reaction norms. Further analysis suggested that clutch size was affected by heterozygosity across the genome, but egg size had more complex relationships, with evidence favouring the influence of multiple loci. Given the apparent lack of inbreeding and large population size, our results suggest associative overdominance as the likely mechanism for the impact of heterozygosity, but also created a puzzle about the process producing associations between neutral markers and the genes affecting clutch size or egg size. One possible explanation is a long-term residual effect of the historical bottleneck that occurred when house sparrows were introduced into North America. The existence of heterozygosity-fitness correlations in a population with considerable phenotypic plasticity and little inbreeding implies that the effects of heterozygosity may be more significant than previously thought.

Individual variation in parental care reaction norms: integration of personality and plasticity.

Personality (consistent differences between individuals in behavior) and plasticity (changes within individuals in behavior) are often viewed as separate and opposing phenomena. We tested this idea by analyzing parental care reaction norms in a bird that exhibits biparental care. Personality in provisioning behavior existed (r(ic) = 0.11) and persisted despite being reduced after accounting for individual differences in environment. Plasticity was also evident and differed between the sexes. Male visit rate was associated with changes in brood size and time of day, but female visit rate was associated with changes in nestling age and date. In both sexes changes in visit rate were positively correlated with changes in their partner's visit rate. Both sexes also exhibited multidimensional reaction norms; interaction terms revealed that within-individual visit rates increased more steeply with brood size when nestlings were older, and the effect of the partner's visit rate was sensitive to variation in date, precipitation, and the focal bird's age. Individuals also varied in how they responded (reaction norm slope) to changes in nestling age and partner visits. Moreover, parental personality was interdependent with individual plasticity in several ways. Individuals of both sexes with a high visit rate also responded more positively to changes in nestling age, and males also showed this pattern with changes in partner visit rate. Explicit use of the behavioral reaction norm integrated personality and plasticity, revealed that these are not opposing concepts, and stimulated new hypotheses about sexual conflict over care and provisioning as a life-history trait.

Heterozygosity and extra-pair paternity: biased tests result from the use of shared markers.

Recent studies of extra-pair paternity have found support for the idea that heterozygous males have an advantage in siring offspring. Most studies use DNA microsatellite loci to determine paternity and then use the same loci to estimate individual heterozygosity. However, because the likelihood of detecting extra-pair offspring depends on the combinations of parental alleles, it is possible that biases arise from particular allele combinations. This might produce false support for the influence of heterozygosity on mating behaviour. We used a simulation model to assess how large this bias might be. We found two sources of bias. First, we found a bias in the null hypothesis of a simple statistical test commonly used to test several predictions of the heterozygosity hypothesis. The use of randomization tests could eliminate this bias. Second, we found that using the same loci for both paternity and heterozygosity can cause an increase in results supporting the heterozygosity hypothesis when no effect of heterozygosity actually exists. This bias is reduced through the use of more markers with higher levels of polymorphism and heterozygosity, but can be eliminated entirely by using a separate set of markers to determine paternity and assess heterozygosity. The two sources of bias reduce evidence favouring the heterozygosity hypothesis, but do not negate all of the studies that support it. We suggest that further studies of heterozygosity and extra-pair paternity are important and likely to be informative, but our recommendations should be incorporated by researchers to improve the reliability of their conclusions.